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SUMMARY 

The laminar hoimdary layer on a circular cone at angle of attack 
to a supersonic stream is analyzed. A solution for the outer nonviscous 
flow, as obtained hy a method of perturbation in angle of attaok, is 
available and is presumed to govern the development of ^the boundary 
layer, subject to certain restrictions and corrections! 

In the limit of vanishing cone angle, the boundary- layer equations 
reduce to those for the flow about a yawed infinite cylinder in com- 
pressible flow, but may not be solved in the simple manner appropriate 
to incompressible flow. 

For small angle of attack, the boundary- layer quantities are con- 
sidered to differ only slightly from the corresponding quantities in 
the known flow at zero angle of attack. ' The system of equations result- 
ing from this assumption are solved to yield the effect of small angle 
of attack on the velocity profiles, skin friction, boundary- layer thick- 
ness, and the forces and moment on the cone, or, more precisely, the 
rates of change of these quantities with angle of attack, evaluated at 
zero angle of attack. 

All boundary-layer quantities show parabolic similarity in merid- 
ional planes . Circumferentially, the boundary- layer thickness associated 
with the profile of the circumferential velocity component is constant, 
and the profile itself exhibits no tendency to separate; -vihereas the 
profile of the meridional component undergoes changes caused by angle 
of attack which show the expected effects of secondary flow in the 
boundary layer. 


INTRODUCTION 

The flow around bpdies of revolution at positive angles of attack 
is characterized by a draining of the boundary layer from the underside 
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to the top, causing the formation of "lohes" of low-energy air near the 
top of the hody. For sufficiently large angles of attack, separation of 
the cross- flow may result in a transversely shed Karmian voirtex street. 
Under these circumstances, large viscous contidhutions to li f t, drag, 
and pitching moment on the hody may be expected. Of course, viscous 
lift, drag, and moment are expected at any angle of attack, whether or 
not cross-flow separation occurs . 

The present investigation, conducted at the HACA Lewis laboratory, 
is concerned with a special case , namely , the laminar bou n dary layer 
"vdilch develops on a id^t circular cone at a small angle of attack to 
the incident supersonic flow. Solutions of the laminar boundary- layer 
equations are sought in order to describe the viscous effects due to 
angle of attack j and particularly those contributing to lift, drag, and 
pitching momeht. 

As is . customary, it will be assumed that the nonviscous solution 
for the supersonic flow about a cone at angle of attack is valid, except 
in a thin viscous layer enveloping the cone surface . The nonviscous flow 
is conical in the Busemann sense j that is, physical quantities are con- 
stant along rays from the apex. According to the analysis of Stone 
(reference l)', the flow quantities may be considered to be the sums of 
the corresponding quantities of zero angle of attack and small increments 
linear in angle of attack. This theory forms the basis for the extensive 
tables prepared at Massachusetts Institute of Technology imder the super- 
vision of Kopal (references 2 and 3) .- The tables of references 2 and 3 
will be used in the present report. In reference 4, Ferri criticizes 
the Stone-Kopal treatment of the entropy gradients near the body and 
proposes a correction to the tables of reference 3. The applicability 
of this correction will be discussed in a subsequent section. 

The boimdary layer itself is to be treated in the manner of refer- 
ende 5, \jherein it is shown that "conical" bodies in supersonic flow 
have boundary layers I'jhich develop parabolical ly along generators of the 
surface. This boundary- layer development is- a consequence of the fact 
that the outer nonviscous flows over such bodies have vanishing velocity 
and pressure gradients along rays from the apex. 

The special case of a right circular cone at zero angle of attack 
has been solved by Hantsche and Wendt in reference 6 and by Mangier in 
an unavailable report. These solutions give a boundary layer which 
thickens parabollcally along cone generators, and provide a formula for 
viscous drag which is valid also for small angle of attack, as will be 
shown subsequently. 

A comparison is made herein between the present analysis and the 
hypothesis, which has been discussed by Allen (reference 7) and Van Dyke 
(reference 8), that the viscous cross flow associated with the fll^t 
of bodies of revolution at angle of attack be regarded as essentially 
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that -Which -would occur in the case of an infinite circ\ilar cylinder at 
yaw -id-th respect to the incident stream, a flow which has heen analyzed 
hy Prandtl, hy R. T. Jones, and hy Sears (reference 9). The attractive- 
ness of this hypothesis lies' in the fact that, for the yawed infinite 
cylinder in incompressihle flow, the cross flow may he obtained from 
-two-dimensional theory without regard to the spamd.se coirponent. In 
reference 5, however, it is pointed out that the chord-wise velocity com- 
ponent does depend to some degree on the spanwise component in the case 
of conrpresslhle flow about a yawed infinite cylinder. Whether or not 
this consideration is of practical importance has not yet been ascertained. 


R03WISC0US FLOWJffiOUT CORE AT ANGLE OF ATTACK 


In order to -breat the boundary layer on a cone at angle of attack, 
it is necessary first to describe the non-viscous flow about the body, 
particularly the velocity components and properties of state -which would 
be predicted at the surface. Figure 1 shows the coordinate system and 
other notation to be used. The distance x is measured along generators 
of the conej y, normal to the cone surface; and <p, circumferentially 
around the cone. The velocity Qomponents u, v, and -w are in the 
X-, y-, and <I>-dlrection, respectively. The cone has a semivertex 
angle of 9 and is at an angle of attack a. , The subscript 0 denotes 
conditions in the xmdlst-urbed stream, and subscript 1 denotes eval- 
uation of the non-vLscous flow at the cone surface, or, alternatively, 
at the outer edge of the boundary layer. The symbol notation used herein 
appears in appendix A. 


Since -the present report is concerned -wi-th small angles of attack, 
and correspondingly Ftmal 1 departures of the entire flow from that occur- 
ring at zero angle of attack, it will be convenient to refer all physical 
quantities to this basic non-vlscous flow. Thus, denoting by bars the 
properties at zero angle of attack, evaluated at the cone surface, the 
foUo-^d-ng symbols on the left -will hereinafter be identified -with 
the dimensionless groups on the right: 


u, w 

P 

T 


u/u, w/u 
p/p 


'u2 


2CpT/ 


~ p/pu2 J 


( 1 ) 


Thus, for example, T ~ 2CpT/u2 
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.On the s\irface of the cone, the theory for small angle of attack 
used in reference 3 yields results "which may he "written in the follo"viing 
form: 

Uj = 1 - a. cos 


= a A 2 sin <P 
= 1 + a A 3 cos (p 


> 


= 1 + a A 4 cos q> j 


( 2 ) 


"Where the A’s correspond as foHo"ws to the quantities (e"val"uated at 
the cone surface) tabulated in reference 3 and do not themselves depend 
on a: 


Ai ~ -x/u 

A2 ** -z/u - 2x/(u sin S) 
A3 - "h/p 

A4 ~ ?/p 


( 3 ) 


The barred (a = O) quantities are tabulated in reference 2. The 
quantities on the right side of eqxiations (3) are in the symbol notation 
of reference 3. 


The relations (2) pro"vlde the boundary conditions to be imposed on 
the "vdscous flo"w at the outer edge of the boundary layer. 


If the definitions used in reference 3 were to be followed strictly, 
then A 2 ~ z/u. The corresponding relation in set (3) is written dif- 
ferently in order to compensate for a systematic error in the tables. 

This error, confirmed by Professor Kopal in a letter to the Lewis 
laboratory, consists in the use of the wrong sign for the quantity C 
appearing in equation (34) of reference 3. Only the values of z (and 
hence Ag) are affected. Since C refers to the deviation of entropy 
in the flow due to the cone thickness, the correction described is of 
little consequence for small cone vertex angles, but is substantial for 
the larger vertex angles. 

In reference 4, it is pointed out that the relations (2) imply that 
"the change in entropy at the cone surface due to angle of attack varies 
circumferentially around the cone as cos <P, "whereas physical reasoning 
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indicates that at the cone surface the entropy should, he constant except 
for a singularity at <p = :t. This difficulty arises because the method, 
of linearization in angle of attack used, in reference 3 is improper near 
the cone surface . In reference 4 it is concluded, that the tables of 
reference 3 are correct, except in a "vortical layer" of an angular 
thickness of order a? (infinitesimally thin in the linear approximation) 
across "which there are a sudden change in entropy of order a and 
corresponding changes in other flow quantities, such as velocity and 
temperature, but not pressure. The "vortical layer" is a transition 
region in which the fona of the entropy variation must change from the 
type cos <p away from the surfeice to constant at the cone surface. 

The corrections to the tables of reference 3 required to take 
account of this effect are not used in the present report; that is, 
the boundary layer will be considered as governed by the flow external 
to the layer alluded to in reference 4. The following argument is pre- 
sented to Jvistify this point of view: The reason for considering the 

entropy constant at the surface is, briefly, that the fluid particles 
nearest the surface must a1 1 have passed through the conical bow shock wave 
In the vicinity of <p = 0 and, once past the shock, will conserve their 
entropy at the value appropriate to the strength of the shock at (p = 0. 
Thus, the sharp normal gradients near the surface (discussed in refer- 
ence 4) may be considered to exist by reason of the conservation of 
entropy along the streamlines nearest the wall. If the boundary layer is 
much thicker than the "vortical layer" (a^), the boundary layer will 
immediately entrain most of the fluid in the vicinity of q> = 0, which ■ 
would otherwise proceed along the surface to form the base of the vor- 
tical layer. The properties of the entrained particles will then dif- 
fuse through the boundary layer by the action of visc’osity, entropy will 
no longer be conserved next to the wall, and the constraint giving rise 
to the vortical layer will be removed. 

Thus, it is reasonable to suppose that viscous diffusion in the 
boundary layer will prevent the formation of the vortical layer described 
in reference 4, and that the boundary layer will be governed at its 
outer edge by the flow described in reference 3 (equations (2)), provided 
that the boundary layer is much thicker than would be the vortical layer. 

If the reverse were true, the boimdary layer would certainly be governed 
by conditions more nearly appropriate to the base of the vortical layer. 

The criterion for entrainment and diffusion of the vortlceil layer would 
be 

8 »xa?. 

where 8 is the boundary- layer thickness. As a rough approximation 
(reference 6), 

1 

8 » (x/r)2 
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■wiiere R is Reynolds number per unit length.. Tb.\is^ the criterion is 

R^« a-4 , ' 

Since the present analysis -will conce3?n only small angles of attack, and, 
strictly speaking, Td.ll he precise only in the limit a-*-0, consistency 

requires that the inequality Rx«<x"'^ he presumed applicable. Of 
course, since the vortical layer is eissociated yri-iii entropy variation, 
this question is of less Importance for smaller cone vertex angles. 


LAMTMR EOUHDARY LAYER ON CONE AT ANGLE OF ATTACK 
Differential ^tmdary-Layer Equations 
and Boundary Conditions 


In reference 5 it is shown that the differential eq^uatlons for the 
conpressihle laminar boundary layer on a conical body in supersonic 
flow may be written: 


[f + 1 g + I gj] fxx - I 8xfX£ + I (8^)" + 2f^XX = 0 (4a) 


[5 I ® + I 8 J Sxx - I Sx«x£ - I hfx - I + 2SXXX = 0 


T + (f^)^ + = T]^ + + -w^2 


P=^i^pT 


(4b) 

(4c) 

(4d) 


Equations (4a) and (4b) are momentum equations, equation (4c) is an 
energy balance, and equation (4d) is the equation of state. The func- 
tions f(X,^) and g(X, ^) axe related to the two-component vector 
potential discTissed in reference 5 and axe defined according to the 
relations 


u = 
w = 



(5) 


and in a manner such as to satisfy the continuity equation identically. 
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These equations embody the following physical assumptions, in 
addition to those arising from the concept of a thin houndary layer 
(for example, that pressure does not vary across. the houndary layer): 


(,a) Prandtl number of. 1, and constant specific heats- 

(h) Uo heat transfer throng the body surface 

(c) Validity of the temperature-viscosity relation 

H = c ? 

p T 


(6a) 


where -C is defined as follows, in order to match relation (6a) to the 
Sutherland formula at the cone surface (denoted by subscript w) : 

1 

C = IfS/T (St) 

V T / t^/t + s/t 

the quantity S being taken equal to (216° R) 2Cp/u2. This temperature- 

viscosity relation as applied to flow over a flat plate is discussed by 
Chapman and Rubesin in reference 10. Since the case of Prandtl number 
of 1 and no heat trsCnsfer is under consideration, T^ may be taken equal 

to the stream stagnation temper atvire (equation’ (4c) ) . 

The coordinates x and y (fig. l) have been made dimensionless 

by referring them to the length C ==. Thus, the normalizing relation 

pu 


pux puy 

x,y - - — , ^ 

Cp Cp 

should be added to the conventions (l) previously established. 


( 7 ) 


The coordinate X has been foimed as follows: 



P dy 


1 

2 


( 8 ) 


Equations (4) and definition (8) show that similarity of the Blasius type 
exists in meridional planes. Thus, for example, boundary- layer thickness 
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is proportional to xV^ and may "be said to increase paratolicaHy 
vith X. This fact is a consequence of the vanishing pressure gradient 
along generators of a conical hody in supersonic flow . 


The variable ^ is a dimensionless coordinate -which, in the case 
of a circular cone, may he replaced hy <p sin 0, or by (p0 where 
0 = sin 9. Equations (4) may thus be rewritten: 




+ 1 P*(<P) g 
30 p ® 





XXX = 0 


(9a) 


S + M S?] ®XX “ ^ - I Sx^x - is + ^Sxxx - 0 

(9b) 

T + (f\)^ + (Sx)^ = T^ + U]^2 + ^^2 (9c) 


P = 



(9d) 


The boundary conditions on the functions f(X^q>) and g(x^q>) 
are; At the outer edge of the bomdary layer, the u and w velocity 
components should take on the corresponding non-viscous values 




(10a) 


gx(»,<p) = Wj^((p) 


(10b) 


and at -fche cone surface, the u, v, and w velocities should vanish, 

fx(0^«P) = gx(0,<P) = 0 (10c) 

f(0,<p) = g(0,<p) = 0 (lOd) 

Equation (lOd) is eq-ulvalent to the requirement that the normal velocity 
component v vanish at the cone surface. This boundary condition is 
discussed at some length in reference 5. 


Simplification of Equations in Special Limiting Cases 

The problem of finding solutions of equations (9) for f and g 
subject to boundary conditions (lO) is rather difficult. Accordingly, 
it is desirable to investigate the manner in -\)hich the equations 
simplify in certain special circumstances. 
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Small cone angle, large angle of attack . - SymboJjLcal ly, tM.8 case 
may "be' descrilDed as follows: 

®«1 

a 

(X of order 1 

These -assun 5 )tlons would "be expected to lead to the equations for the 
boundary layer on a yawed infinite cylinder (the conpressihle analog 
of the flow treated in reference 9) . For the purpose of this dis- 
cussion, it is convenient to define new variables: 


A = 


= X _ 





p <3y 


(0x) 


_1 

2 


(11) 


F H f/^^ 

G = g/yie 

Equations (9) and (lO) become 

[sep + PiM a + 2 gJ Fjii - 2GjiTi, + 26 


2F 


= 0 (12a) 


*3 


30F + 


P*(<P) 


G + 2G(^ G aa - 2G^Gj^(p - 20(^F^ - 2 
5a(«^<P) = Ui(<p) 

Ga(«^<p) = Wj^(q>) 

Fa(0,<p) = Gj^(0,q>) = F(0,(p) = G(0,<p) = 0 


P'(<p) 


+ AAA ® 


(12b) 


(13) 


For flow about a cone at angle of attack, F and its derivatives will be 
presumed to be of unit orders G and its derivatives, and the derivatives 
of P 3 _, of order a (order 1, in this case). When terms of order 0 
are neglected, equations (l2) may be written; 


^A^A(p " ^AA - ^2p^ ®"aA + 


^AAA 


(14a) 


- ^<p^AA = - (1 - % GGaa) + GaAA (14-b) 
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Equations (li) are also the equations for the compressihle visco\is 
flow about a yawed infinite cylinder and, for the inconipressihle case, 
are precisely equivalent to those treated in reference 9. Inspection 
of definition (u) shows that similarity exists in terms of the reciprocal 
sq.uare root of the local radius (6x) of the cone cross section'. For 
an infinite cylinder at yaw, this radius would he a constant. For a 
sufficiently slender cone at a sufficiently large angle of attack 
(e/a«l), the hypothesis discussed in reference 7 is therefore essen- 
tially correct j namely, that at any point on the cone, the boundary 
layer is that "vdiich would appear for the same angle of attack on an 
infinite circular cylinder having a diameter equal to the cross- 
sectional diameter of the cone at the point under consideration. 

Equations (l4a) and (l4b) are coupled through the density p, -which 
appears in the pressure gradient term, and therefore G cannot be 
obtained independently of F as in the incompressible case. Compress- 
i-bility enters in this way according to the stream Mach nuniber. In 
reference 8, it is suggested that the important Mach number is that 
based on the cross velocity. 

This limiting case retains the general feature of parabolic gro-wth 
of the boundary layer along cone generators despite the fact that, in 
the limit, the configuration is -that of a cylinder of infinite length. 

Small angle of attack, large cone angle . - This situation, may be 
specified as 


to 

CO 

c\a 

CV] 


6 of order 1 

presumably yielding, in -the limit, the equations for the boundary layer 
on a cone at zero angle of attadk to a supersonic stream. Equations (2) 
show that for small angles of attack, u-j^(cp) » 1 and w-j_(<p) and p'(<P) 

are of order a. Therefore, from equations (lO) , f (X,<p) will be of 
order 1 and g(k,<p) -^d-ll be of order a. With terms of order a in 
equations (9a), (9b), and (lO) neglected, equations (9a) and (lO) become 

+ 2fxXX = 0 (15a) 

f^(0,cp) = f(0,q>) = 0 (15c) 

Equation (9b) drops out because 1 its terms are of order a or smaller. 
In -view of the foim of equations (15), f is a function only of X, and 
equation (l5a) may be interpreted as an ordinary differential equation. 
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Eq.uations (15) govern the houndary layer on a cone at zero angle of. 
attack to a supersonic stream and are formally Identical with the equa- 
tions for plane flow over a flat plate, the difference "between the cone 
and flat plate cases appearing only throu^ a slightly different definition 
of independent variable. (See, for example, reference 6 or Mangier 's 
work.) 


M 

CO 

CD 

CD 


PERTUEBATION ANALYSIS FOR SMALL ANGLE OF ATTACK 

In the remainder of this report, the laminar houndary layer on a 
cone at angle of attank to a supersonic stream will he analyzed as a 
perturbation of the flow at zero angle of attack (equation (15)); that 
is, under the conditions 


e 


<1 


0 of order 1 


Formulation of Perturbation Equations 

With the assumption that the entire flow may be represented as the 
sum of the basic flow fo(^) satisfying equations (l5) and correction 
terms proportional to angle of attack, f (X,cp) and g(X,(p) may be 
written 

f = ^0^^^ - a A^ cos <p f^(X) + . . . (I6a) 

g = a Ag sin <p g 3 _(X) + . . . (l6b) 


The quantities fg, f]_, gj^, and A 2 are presumed of order 1. The 

forms of the correction terms are chosen to be consistent with the 
relations (2) to be imposed as boundaay conditions on the boundary layer 
at its outer edge. 


The pressure gradient terms appearing in equations (9a) and (9b) 
may be expressed in perturbation form with the aid of equations (2), 
(9c), and (9d) and definitions (16): From equations (2), since pressinre 

is presumed not to vary across the boundary layer. 


P 


-a Aj sin (p 


(17 a) 


From equations (2), (9c), and (9d), 
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1 _ r -1 T _ r -1 ^1 ^1^ ^1^ ~ “ (sx)^ 

p (l + a A 3 cos (p) 


2 r p 


2r 


(1Tb) 


Because differs from T only by a quantity of order equa- 

tions (I 6 ) may be substituted into equation (l7b) to yield, for 
small a. 


P 

P 


r -1 

2r 


+ 1 - (fo’) 

tdiich may be multiplied into equation (l7a) to give 

i2 


S;M = -a JW sin » ni 
P 6 2f 


T + 1 - (fo>) 


(17c) 


Substituting equations (I 6 ) and (l7c) into (9a) and (9b) and 
equating the sum of terms of unit order in equation (9a) to zero give 


Vo" 2fo"'= 0 


(18) 


The sum of terms of order a in eqtiation (9a) is next equated to zero, 
yielding 

* “i”'" ° 


The largest terms in eqxiation (9b) are of order a; the sum of these 
terms is equated to zero: 


2^oSi" - 2^0 'Si' 6gi' 


= -2 


^5 r-1 

9A2 2r 



^ - (f, 



( 20 ) 


Substituting equations (I 6 ) and (2) into boundary conditions (lO) 
and equating terms of like order in a provide the boundary conditions 
to ■which the differential equations (18), (19), and (20) are subject: 


fo«(«) = 1 

fo‘(o) 

= = 0 

(21) 

f^>(«) = 1 

fi'(o) 

= fi(o) = 0 

(22) 

H 

II 

T 

gl’(o) 

II 

o' 

II 

0 

(23) 


2286 



qR23 


KAjOA TU 2521 


13 


In view of iDonndary conditions (2l) and (23), the relation 


T = 

QAg 2 


1 


(24a) 


must hold in order that eq.uatlon (20) he satisfied at X = “>. 
of relations (l). 



Because 


(24h) 


where M is the Mach numher at the surface of the cone when at zero 
angle of attack. In terms of guantlties tabulated in reference 3, 


M^ = 


u^ 


r-1 i-u^ 


(24c) 


With the aid of equations (3), (24h), and (24c), it may he verified from 
reference 3 that equation (24a) is indeed correct. Equation (20) may 
therefore he written more simply as follows; 


3fQg^" - Sf^'g^’ + 6g^«“ = -2 \1 + M^ 


- 


.^2 


(25) 


To summarize, the basic flow at zero angle of attack over the cone 
is determined by the solution of the boundary-value problem 


Vo" + 2f^’” = 0 


fQ«(») = 1 f^*(0) = f^{0) = 0 


( 18 ) 

( 21 ) 


The merldlonEa velocity component u is given by the function fQ’(X), 
according to equation (5) . When fo(X) is known and M is specified, 
the circumferential velocity component due to a small angle of attack 
may be determined from equations (5) and (l6b) and from the solution 
of the problek ' 


3fQg^" - 2fQ«^« + 6g^”' = -2 <1 + ^ 


.^2 


1 


(25) 


g^«(«) = 1 


g2*(o) — gj_(o) — 0 


(23) 


Then, presinnlng fQ(X) and gjL(l-) to have been determined and the 
quantity A^/QA-^ to be specified, the Increment in meridloneil velocity 



14 


MCA TPr 2521 


component due to a small angle of attack may "be found from eq.uatlons (5) 
and (l6a) and from tlie solution of the problem 

3f0fl" + 3fo"fi + 6fi“' = 2|^g^o"gi (19) 

f^'(««>) = 1 f^«(0) = f^(0) = 0 (22) 

The perturbation procediire used to obtain these differential equa- 
tions is, in general, subject to an important limitation -which should 
be examined. Since pH velocity components vanish at the cone surface, 
it is perhaps not ob-vious that it is proper, near the surface, to make 
the neglectlons in equations (9) necessary to arrive at equations (19) 
and (25) . In particular, consider the term ff\x appearing in equa- 
tion (’9a) . Substitution of equation (l6a) yields 

ffX\ = fo^o" - cos <p (f2_fo" cos^ q> fif^" + . . . 


(26a) 

In deri-vdng equation (19), the term involving has been neglected 

in comparison -^Tith the term involving a, despite the fact that both 
terms vanish as X-K). The problem of proper neglection near the sur- 
face becomes determinate if the last two terms of expression (26a) 
are written as follows: 


( f f ” 

-a A]_ cos (p (f^fo" + il - a cos (p 

V f.f^" + f, 


t .p ^ It 


Thus, -fche question is -whether or not, for a« 1, the folio-wing Inequality 
holds as X-»-0: 


a A^ cos <p 


flfl" 


f f " + f f " 
10 0 1 


«1 


(26b) 


or, equivalently, whether or not 


Vi" 


(26c) 


is always of order 1 or smaller. 
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Boundary conditions (2l) and (22) indicate that Tq and may 

he represented for small ^ hy the following power series; 

^0 = ^ + FT • 

^ ^ + . . . 

If fQ"(o) ^ 0 (as in the present problem), then, in the limit \->0, 

quantity (26c) is equal to i f 2 ^”( 0 )/fQ"( 0 ), which is presumably no 

2 

larger than nalt order. All the neglections made in deriving eqna- 
tions (19) and (25) may be justified by this sort of argument. 

If, however, fQ"(o) were zero and f2^”(0) were not zero, expres- 
sion (26c) would become infinite as X-*-0, and Inequality (26b) would not 
hold. The quantities fQ"(0) and f-j^”(o) are related to the slopes 

at the wall of the basic and incremental velocity profiles (equa- 
tion (5)), and hence are proportional to the corresponding contributions 
to viscous shear stress at the wall. Thus, in general, a pertu’-bation 
analysis of a bound^y layer will be improper near a point -vdiere the 
shear stress of the basic flow vanishes, that is, near a separation 
point . 

Solution of Perturbation Equations 

The boundary- value problem for Tq (equations (18) and (2l)) arises 

in the case of boundary- layer flow on a flat plate and has been studied 
thoroiighly in the literatirre. For conpleteness, the values of fQ and 

its derivatives, as determined by calculations made at the Lewis 
laboratory, are presented in table I. The results agree precisely with 
those presented in reference 9. 

The functions f]^(^) and g]^(X) may be expressed as linear com- 
binations of functions not depending on the parameters M or A^/OAj^. 

The following representation was proposed by L. Richard Turner of the 
Lewis laboratory; 

gj^(\) = fQ(X) +^1 + M^)h^(X) 


(27a) 
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(X)= (2 Mf (X) + 1 - 2 MhgCX) + 2(1 + Ili li 3 (X) 

\3 9AJ V 3 6aJ \ 2 /0Ai 


(27b) 

The coefficients in equations (27) may be determined from tbe tables of 
reference 3. Substituting equations (27) into equations (19), (22), 
(23), and (25) and making use of equations (18) and (2l) yield tbe 
following differential equations and boimdary conditions for tbe func- 
tions and bg, wblcb can be solved once and for all, witb no 

regard for tbe physical parameters of tbe problem: 

3fQbi" - 2 fQ«bi* + 6 b 3 ^«'« = - 2 |l - (fQ*)^ (28a) 

b^«(») = b^'(o) = b^(0) = 0 

3fQb2" + 3fo"b2 + 6b2"‘ = 0 

b2«(»)=l b2«(0) = b2(0) = 0 

3f b " + 3f ”b + 6 b • " = f "b 
0 3 0 3 3 0 1 

b3>(«) = b3«(0) = b3(0) = 0 

In appendix B, tbe computation of bj_, b 2 , and bj is described. Tbe 
results appear in table I. 


(28b) 

(29a) 

(29b) 

(30a) 

(30b) 


Presentation and Discussion of Results 

Tbe various important boundary- layer quantities may be found from 
tbe profile functions given in equations (27). Tbe numerical examples 
presented in figures 3 to 15 bave been computed witb tbe aid of tbe 
tables of references 2 and 3. 

Tbe boundary layer on a cone at angle of attack is cbaracterized 
by tbe phenomenon of "secondary flow," -wblcb has been described by 
various authors. A discussion of secondary flow as it occurs on a cone 
at angle of attack is pro-vlded in tbe folio-wing paragraphs in order to 
facilitate interpretation of tbe results of tbe present analysis. 


«?286 
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Since the Inviscid flow about the cone at angle of attack is 
"conical," the pressure gradient at the outer edge of the houndary layer 
is entirely circumferential (equations (2)) . !Ehe streamlines in the 
Inviscid flow jiist outside the hotaidary layer (shown schematically eis 
solid lines in fig. 2) will incline toward the direction of the pres- 
sure gradient. Since the presstire tends to he constant across the 
houndary layer, the fluid in the boundary layer is subject to the same 
pressure gradient as is the outer flow, hut has less Indrtia with which 
to resist its effect, and thus it tends to follow the direction of the 
pressure gradient more closely than does the outer flow. 

An angular divergence therefore exists between outer streamlines 
and houndaary-layer streamlines; the divergence is greatest when the 
limiting streamlines at the base of the boundary layer are considered. 
These limiting streamlines are sho\m schematically as dashed lines in 
figure 2. The divergence 6 Is defined as the local angle Included 
between the outer and limiting streamlines. 

This divergence between inner and outer streamlines ("secondary 
flow") clearly leads to a draining avray of low-energy air from the 
hl^-pressiire side of the body and a corresponding concentration of 
low-energy air on the low-pressure side. Accordingly, if the cone is 
at positive angle of attack, the skin friction would be expected to be 
lower at the top of the cone than at the bottom, and the boundary- layer 
thickness would have a maximum at the top and a mlnlnuTm at the bottom. 

Before proceeding firrther, it is perhaps wen to emphasize that the 
present analysis is valid only in the limit of vanishing angle of attack. 
That is , terms multiplied by a represent rates of change irtth angle of 
attack, evaluated at zero angle of attack. All results win be presented 
in this form. Whether or not the absolute change for a sman finite 
angle of attack can be obtained from this theory depends on the relative 
mEignItude of the effect so computed, rather than on the size of a, nnt^ 
depends further on the (unknown) second and higher derivatives of the 
quantity with respect to a. 

Velocity profiles and streamline divergence . - From equations (5) 
and (16) , 

w = gx = a Ag sin (p g 3 _' 

u = = f^’ - a A^ cos,(p f^* 

These expressions may be divided by the values of 
in equations (2) to yield, for a«l^ 


(31) 


W]_ and Uq^ given 



(32a) 
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^ = fo' + a COB (p (fQ* - f^') 

■Ui 

or, substituting equations (27), 


(32t) 


(33a) 


— = + a xi^ cos <p 

U]^ ^ J- 


1 - ^ 


2 ^2) 


3 0A-, 


(f^« - h_') -21+ 


--1 zA A2 


M 1 — ^ lu' 
2 /©A^ ^ 


(33b) 

The outer streamline is inclined to the x-direction hy an angle 




-1 W, 

tan"-^ = tan' 


■If ^ ^ 

^1 - a Aq^ cos <py 


a A 2 sin (p 


which is obtained by ^xse of equations (2). The limiting streamline 
is inclined to the x-direction by an angle 


tan 


-1 


(lim Z 1 » a Ap sin <p 
Vx-K) uy ^ ± q '' 


Si" jo) 
( 0 ) 


which is obtained "(Tith the aid of equations (3l) and L 'Hospital’s 
rule. The divergence e is given by the difference between these two 
angles : 


e » a Ag sin <p 
or, by use of equation (27a) and table I, 


g "(0) 

_± - 1 


^q"{ 0) 


e « a(l.29l) Ag ll + ^ M 1 sin q 


-i\ 


I 


(34) 


This streamline deflection due to angle of attack is indicated by the 
quantity 
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— L_ = (1.291) A2 fl + ^ 

sin (p \Sa/cfcO \ 2 J 

and is slao'vjn in figure 3 as a function of stream Mach numher Mq for 
cone semivertex angles 0 of 10°^ 20°, and 30°. 

The dLstrihution throng the houndary layer of the circumferential 
velocity -w/w-j^ is shown in figure 4 for cone semivertex angles 0 of 

10° and 30°, and for stream Mach numhers Mq of 2 and 4. From equa- 
tion (3l) and figure 4, it is evident that no "cross-flow separation" 
takes place, because the profiles of figure 4 (tdiich are not of the 
separation type) appear at all circumferential locations, modified by 
the value of sin <P, ivhich affects only the scale of the w-prof iles . 

Of course, cross-flow separation TOuld occur for large angles of 
attabk because of the effect of an adverse circumferential pressure 
gradient over part of the cone sinrface. For very small angles of 
attack, separation of the cross flow would hardly be expected since the 
circumferential pressiare gradient is always favorable (equations (2)).. 

/ 

The increment in meridional velocity profile due to angle of attack 
is indicated by the quantity 

°°® 'P feUo^O = \ "^l' ^ 

which is sho\m in figure 5 for 0 = 10° and 30°, and for -Mq = 2 and 
4. For positive angle of attack, the meridional profile is correspond- 
ingly steeper on the underside and less steep on the top. This result 
agrees with the queilltative argument that low-energy air should drain 
aimy from beneath the body and concentr^-te on the top. The order of 
magnitude of this effect is shorn in figure 6, in -which the n./u^^ pro- 
files in the plane <p = 0,rt for a = 2.5° are compared with the pro- 
files for a = 0, when 0 - 10° and Mq = 4. 

The relative thickening of the boundary layer at the top of the 
cone «nri thinning at the bottom might be expected to have an adverse 
effect on laminar stability at the top and the reverse effect on the 
bottom. 

Skin friction . - The clrcvimferential and meridional components of 
the viscous shear stress at the cone sirrface may be written as 
coefficients; 
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•where "the quantities making tip "the ri^t sides are in dimensional fonn. 
Application of equations (l) and (4) to (8) yields 



(34a) 


(34h) 


•where 16 a Reynolds number based on distance from the cone apex 

and on fluid quantities e^vnluated in the non'viscous flow at the surface 
of the cone at zero angle of attack: 


Rx 


pu X 


Introduction of equations (2) and (16) yields 




A 2 sin <p g 3 ^"( 0 ) 



a 


A^ cos q> 


/ ^l”(0) 

VV(o) 



or, using equations (27) and table I, 



a 




sin 


<p^l. 


522 + 0.858 


--1 -A 


(35a) 
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0.664 


1 - a cos (p| 1.502 - 0-525 


0An 


0.190 





9A^ 2Aj_ 


(35b) 


The meridional skin friction at zero angle of attack is obtained 
by Betting a = 0 in equation (35b), yielding a result equivalent to 
the results of Hantsche and Wendt (reference 6) and Mangier. 

In figure 7, the circumferential skin friction appears as a func- 
tion of stream Mach number for semi vert ex angles of 10°, 20°, and 30°. 

Figure 8 shows the Increment due to angle of attack of meridional 
skin friction plotted as a function of Mq for the same cone angles. 

The effect of angle of attack on meridional skin friction is of the same 
order of magnitude as the effect on circumferential skin friction, as 
would be eapected in view of the previous discussion of the velocity pro- 
files. For positive angle of attack, the meridional skin fhiction has a 
m a x imu m at <p = 0 and a minimum at q> = jt, which is consistent with the 
previous qualitative -discussion of secondary flow. 

Boundary- layer thickness . - The boundary-layer thicknesses are 
expressed aa the mass-flow defects ( "dlsplaceinent thicknesses") 

* associated with the u and w profiles . Thus 



or, making use of equations (l), (7), and (8) and expressing 6 as a' 
Reynolds number Eg = p u S/p, 
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Substituting equations (5) and (37) into (36), 
1 


5 


no# 


WO 


p/p 2 






2 w. 


dX (38a) 



p/p 2 ^2 


1 - 


u‘- 


u 


u^^y ui 


dX (38b) 


Introducing equations (2) and (33a), performing the indicated integrations, 
and making use of table I cbange eq.uation (38a) to 


= 0.944 + 1.608 (39a) 

It iTOuld be improper to present tbe term of order a resiiltlng from 
integration of eq.uation (38a) since in tbe present approximation tbe 
ratio w/v-j^ is known only to unit order (see eq,uation (33a)). Intro- 
ducing equations (2) and (33b) into equation (38b) and performing tbe 
required integrations (-vritb tbe aid of equations (18), (19), and (25)) 
result in tbe follo-vrlng expression: 
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Rg^/1 -I- = ^1.721 + 2.385 M^j 1 + a cos <p 


2Aj_ A]^ 


0.863 - 


1.721 + 2.385 ^ 


3.576 ^ 


'0.863 + 



The displacement thickness associated with the u-component at zero 
angle of attack (obtained hy setting a = 0 In equation (39h)) Is shown 
In figure 9^ as a function of for various cone angles, and Is 

equivalent to the results of Hantsche and "Wendt (reference 6) and Mangier . 


The w-thlckness given hy equation (39a) Is shown In figure 10. This 
thickness does not depend on <p, and hence, at any value of x. Is con- 
stant around the cone. 


The Increment In u-thlckness due to angle of attack (fig. ll) 
Indicates thinning of the boundary layer on the hlgh-pressiure side of 
the cone and thickening on the low-pressure side. This effect Is qual- 
itatively consistent with the previous discussion of secondary flow In 
the boundary layer of a cone at angle of attach, and Its order of 
magnitude Is consistent with the corresponding results for velocity 
profiles and skin friction. 

Forces and moments . - The Invlscld solution for flow about the cone 
at an angle of attack provides the forces and moments due to the action 
of normal pressiire forces, which may be evaluated using the tables of 
references 2 and 3. In order to obtain the total force and moments, the 
forces and moments due to viscous shear at the surface must be added to 
the corresponding pressure forces and moments. Lift and drag are measured 
normal and parallel, respectively, to the stream direction. Pitching 
moment Is measured about the cone apex, positive In the sense shovm 
In figure 12. All quantities appearing In the following formulas are In 
dimensional form. 
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(a) Viscous drag: Written in coefficient form, the drag due to 

skin friction at the. surface is 



D, 


1 2 . 

1 p u A 

2 0 0 


nx 


rZTt 


1 

A 


— 2 
pu 


PQ^O 


dx 




0^ X dtp 


to 


■VThere D is the viscous drag and. A is the cone surface area. Using 
equation^ (35h) yields 



0.885, 


1 


— - -2 
30 P M 


^ P| 


0 Mo 



(40) 


This expression gives the drag at zero angle of attack. The correction 
due to , angle of attack is of order a?, and hence is neglected in the 
present linear analysis. The drag coefficient (equation (40))^ ■which 
is equivalent to the results obtained hy Hantsche and Wendt (reference 6) 
and Mangier, is shovjn in figure 13 as a function of Mq for various 

values of 0 , 


is 


(h) Viscous lift: The lift due to circumferential skin friction 


nx 








= -a Op +- 
2 ^v A 2 


P ^ 
0 0 


P2jt 


dx 


Of 0x sin(p d(p 
<}> 


0 


■vdaere L is the viscous lift. Using equation (35a), 


V 


C = a 


V 




015 + 0.572 ^ ^ 


2) p_ JEc 
^0 Mq2 k 


(41) 




The lift coefficient is sho'wn in figure 14 as a function of Mq for 
various values of 0 . Except for quite low Reynolds ntunhers, the laminar 
•viscous lift is small conpared with the press'ure lift. The pressure 
lift coefficient is of order a, "^diereas for a Reynolds number based on 
slant hei^t of 300,000, the viscous lift coefficient on a 10° cone at 
= 4 is only of the order of 0.01 a. 
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(c) Viscous pitching moment: By use of equation (35a) ^ the mbment 

due to circumferential skin friction may he found: 



HV 


1 2 . 

5 Vo ^ 



(42) 


■where m^ is the ‘viscous pitching moment about the apex. The moment 
is sho-wn in figure 15 as a function of Mq for several cone angles. 


CONCLUSIONS 

For small angles of a"fctack, the outer nonviscous flo-w about a cone 
at angle of attack to a supersonic stream has pre"vlously been treated 
by a perturbation theory, ■which may be used to describe the outer flow 
governing the developmen-fc of the laminar boundary layer on the cone 
surface, pro'vlded that no accoimt need be taken of the "vortical layer" 
recently dlsciissed by Ferri. 

Gkiis vortical layer is a region in the non'viscous flow field near 
the surface of the cone across which there are large gradients of 
entropy and veloci'ty, and is said to have an angular thickness of the 
order of the square of the angle of attack. If the boundary layer is 
much thicker than this, then presumably -viscous diffusion disperses the 
vortical layer, -vdilch then may be neglected. 

Since the outer non'viscous flow is "conical, " use may be made of 
pre-vlously derived differential eq-uations for a boundary layer governed 
by an outer flow "with conical symmetry. These equations imply a solution 
ha"vlng parabolic similarity in meridional planes . Assuming a Prandtl 
number of 1, no heat transfer through 'the surface, and validity of a 
linear temperature-'vlscosity relation, it is necessary only to solve the 
differential equations for -two components of the vector potential. 

In the limit of vanishing cone angle (finite angle of attack), the 
boundary- layer equations become precisely those governing the compressible 
flow over a yawed Infinite cylinder of diameter equal to the local diameter 
of the cone at the point under consideration. The cross and spanwise 
flows are coupled through the density, however, and it is not proper to 
compute the cross flow by two-dimensional methods, as may be done in the 
incompressible case. This coupling depends on the stream Mach number 
and not on the Mach n-umber of the cross flow. 
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In the limit of vanisMng angle of attack (finite cone angle ) , the 
houndary-layer equations reduce to those governing flow about a cone at 
zero angle of attack. A sufficiently small angle of attack may be con- 
sidered to impose small perturbations on the flow at zero angle of attack. 
If. these perturbations are written in a form appropriate to the stream 
boundary conditions provided by the nonvlscous pertvirbatlon theory, 
ordinary differential equations resiilt. These may then be solved by 
numerical integration to obtain the fluid velocities and state quantities 
in the boundary layer. The type of perturbation method used in this 
report woiILd be Improper in the vicinity of any separation point of the 
basic flow. No such difficvilty arises in the present application. 

I 

According ‘to the perturbation analysis, the boundary-layer thick- 
ness associated ^ri-th the cross (circumferential) velocity is constant 
around the body at a given distance from the apex. The shear associated 
with this cross flow gives rise to lift and moment on the body. No 
separation of the cross flow occurs, presumably as a consequence of the 
favorable clrciomferentiail pressure gradient. 

Secondary flow occurs in the bormdary layer of a cone at angle of 
attack; the slow fluid near the base of the bormdary layer tends to 
incline more toward the direction of the circumferential pressure 
gradient then- does the outer flow. This mechanism causes a concentration 
of low-energy air at the low-pressure side of the cone (top, for positive 
angle of attack) . Therefore, taking into account the change in the pro- 
file of meridional velocity due to angle of attack, the boundary-layer 
thickness associated with this profile is greater at the top and less on 
the bottom of the cone than the value ap;^opriate to the case of zero 
angle of attack. The reverse is true of meridional skin friction. 
Meridional skin friction is Integrated to give the viscous drag of the 
cone. The effect of angle of attack on drag is second order in angle 
of attack and is ignored in the present analysis . 

The changes in the meridional velocity profile due to angle of 
attack suggest that a decrease in laminar stability on the Ibw-pressure 
side and the reverse effect on the high-pressure side migit be ezpected. 

Leiris Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, August 10, 1951 
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APPENDIX A 


A. 


% 


■'nv 


D. 


SYMBOLS 

The following symbols are used in this report: 

related to perturbation in meridional velocity u outside 
boimdary layer 

related to perturbation in circumferential velocity outside 
boundary layer 

related to perturbation in pressure outside boundary layer 

related to perturbation in density outside boundary layer 

constant of proportionality in ten5>erature-viscosity 
relation 

coefficient of viscous drag 
coefficient of viscous lift 
coefficient of viscous pitching moment 
specific heat at constant pressure 
viscous drag 

component of vector potential related to meridional 
velocity u 


^ 0 ^ ^1 


g 




value of f for a = 0 and perturbation of f , 
respectively ■ 

component of vector potential related to circumferential 
velocity 

linear approximation to g 


^1^ ^2'' ^ fimctions composing f^ and g^ 



28 


HACA TU 2521 


M 

m 

V 

P 

R 

T 

u, V, -w 


X, y 


a 


r 

8 



e 


. 0 , 0 


\ 


■viscouB lift 
MacA mmiber 

viscous pitching moment about cone apex 
pressure 

Reynolds number, subscript indicates length upon which based 
temperature 

meridional, normal, and circumferential velocity component, 
respectively 

coordinates taken along cone generators and normal to cone 
surface 

angle of attack 

ratio of specific heats 

boundary-layer thickness 

displacement thickness associated with u-profile 
displacement thickness associated with w-profile 


divergence between outer and limiting streamlines 
dimensionless coordinate («p sin 0) 

semivertex angle of cone and sine of semivertex angle 


dimensionless coordinate 



(j. viscosity coefficient 

p density 

(p angular coordinate measirred circumferentia l 1 y around cone 
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Subscripts : 

V denotes evaluation at cone surface 

0 denotes evaluation in undisturbed stream 

1 denotes evaluation at outer edge of "boundary layer 
Subscript notation for partial differentiation is used. 

Superscripts : 

’ Primes denote ordinary differentiation, 

— Bars (as in u or e) denote evaluation in nonviscous 

flow next to cone surface in. case of zero angle of 
attack. 
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APPENDIX B 

NDMEKECAL DETERMINATION OF PROFILE FUNCTIONS 
By L. Richard Turner 

The four functions fQ(^)j h^(X), hg 

obtained hy n\imerical solution of the houndary-value problems set forth 
in equations (18) and (28) to (30) . 

O 

Differential eq.uatlon (18) is the basic Blasius equation; it has 
been solved to provide the values of the function foC^) req.uired for 

the remaining integrations. Boundary conditions (2l) are specified at 
tuo points X = 0,«» and thus cannot be satisfied directly by a step- 
by-step integration process . This difficulty may be avoided as follows : 
Define 

(l(a) = afQ(X) 



a 


•where a is an arbitrary constant. Substitution into equations (18) 
and (2l) yields 

qq” + 2q”» = 0 (B1) 

q‘(») = a2 q'(0) = q(o) = 0 (B2) 

Differential equation (Bl) may be integrated using the boundary con- 
ditions q'(0) = q(0) = 0 and an eirbltrary value of q"(o). A constant 
value of q'(cr) -will ultimately be reached and may be equated to a^ 
since a^ is arbitrary. The resiiltlng value of a may then be used to 
Invert the transformation leading to equations (Bl) and (B2), thus 
determining the solution fQ(X) . 

Equations (28) to (30) are linear in their dependent variables; 
solutions that satisfy the boundary conditions at X = ® may therefore 
be found by appropriate linear combinations of the complementary and 
particTilar solutions satisfying the given boundary conditions at the 
origin. 

The functions h-j_, h. 2 } and hj were found by computing the 

coirplementary and particulEir integrals of equations (28a) and (30a) 
subject to the follo'Wing boundary conditions: 


(X), and b 2 (X) have been 
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Equation 

Integral 

h(0), h‘(0) 

h"(0) 

(28a' 


Complementary 

0 

1/3 

(28a 


Particular 

0 

1/3 

(30a 


Complementary 

0 

1/3 

(30a, 


Particular 

0 

0 


The values of h"(o) are arhitrary. 

The function h^ Is thus composed of the complementary and 
particular solutions of eq.iiation (28a), linearly comhlned in such a 
way that h^’(®) =0 (boundary condition (28b)). The function h 2 (X) 

was obtained by scaling the con?)lementary integral of equation (30a), 
which satisfies equation (29a), so that h^(«>) = 1 (boundary 

condition (29b)). The function h 5 (X) is composed of the complementary 
and particular integrals of equation (30a), combined linearly so that 
h 2 *(®) = 0 (boundary condition (30b)). 

The four differential equations (l8a) and (28a) to (30a) were 
integrated mnnerically by a method some'what similar to that of Adams 
( see reference H) . Given the complete solution for the function in 
question at six equal 1 y spaced points, the solution may be extended to 
the next point as follows: A sixth-degree polynomial is passed through 

the six known and one unkno^m values of the third derivative. This 
polynomial may be integrated to give the value and first two derivatives 
of the function at the seventh point in terms of the unknown third 
derivative. Substitution of these results into the differential equation 
evaluated at the seventh point yields an algebraic equation for the 
third derivative at that point, •which may be solved, thus extending the 
► solution of the differentieil equation by one step. This process is 
re’peated out to sufficiently large values of X so that -the boundary 
conditions at X = o» may be satisfied to the desired degree of accuracy. 
A step size of 0.01 was used for equation (18a), and a step size of 0.02 
■was used for equations (28a) to (30a). 

The use of this method requires that an initial set of values of 
each function and its first t-\ro derivatives be kno'wn in advance for the 
first six intervals, starting at the origin. These values were obtained 
by use of a method of successive approximation -which is a numerical 
application of the classical Picard iteration process. 

All the numericeil integrations were performed on the IBM Card 
Programmed Electronic Calculator, and carried out to ten decimal places. 
The computations were checked by computing the third and seventh 
differences of the third derivatives to assure freedom from computational 
error and from significant truncation error, respectively. 
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D 


'o' 


1^9 

hi- 

hi" 



hi" 

h3 

^ 3 ' 


0 

0 

0 

0.332056 
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Figure 3. - Streamline divergence acux)8s "boundary layer. 
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(a) SemlTertex angle, 10°. 
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(b) Bemi vertex angle, 30°. 


Tlgure 4. - Profile of clrcumf erenti al velocity component. 
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Figure 8. - Increment in meridional skin friction due to angle of attack. 
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Figure 10. - Displacement thickness of circumferential velocity 
profile due to angle of attack. 



profile due to angle of attack. 
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Elgure 12. - CoAvEirfclons for aerodjuamic forces 
and nsoments on cone. 
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Figure 13. - Coefficient of vIscoub drag. 









